Precise mass-dependent QED contributions to leptonic g-2 at order ex 2 and a 3 
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Improved values for the two- and three-loop mass-dependent QED contributions to the anomalous 
magnetic moments of the electron, muon, and r lepton are presented. The Standard Model pre- 
diction for the electron (g — 2) is compared with its most precise recent measurement, providing a 
value of the fine-structure constant in agreement with a recently published determination. For the 
t lepton, differences with previously published results are found and discussed. An updated value 
of the fine-structure constant is presented in "Note added after publication." 
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I. INTRODUCTION 

The QED part of the anomalous magnetic moment ai — 
(gi — 2)/2 of a charged lepton I = e, // or r arises from 
the subset of Standard Model (SM) diagrams containing 
only leptons and photons. For each of the three leptons 
I, of mass to/, this dimcnsionless quantity can be cast in 
the general form [1] 



mi 



mi 



A 1 + A 2 — + A 2 — +A 3 — — , — - , (1) 
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where mj and mk are the masses of the other two lep- 
tons. The term A\, arising from diagrams containing 
photons and leptons of only one flavor, is mass and fla- 
vor independent. In contrast, the terms A 2 and A3 are 
functions of the indicated mass ratios, and are generated 
by graphs containing also leptons of flavors different from 
I. The contribution of a lepton j to af ED is suppressed 
by (mf/m 2 ) if mj 3>m;, while it contains a logarithmic 
enhancement factor ln(m;/TOj) if mj <C mj. The muon 



contribution to 



is thus power suppressed by a fac- 



tor (m 2 e /m 2 ^) ~ 10~ 5 ; nonetheless, as we will discuss in 
sec. II, this effect is much larger than the tiny uncertainty 
very recently achieved in the measurement of a e [2] . On 
the contrary, the QED parts of a MjT beyond one-loop are 
dominated by the mass-dependent terms. 

The functions Ai (i — 1, 2, 3) can be expanded as power 
series in a/ir and computed order-by-order 

(2) 

Only one diagram is involved in the evaluation of 
the lowest-order (first-order in a, second-order in the 

electric charge) contribution - it provides the famous 

(2) 

result by Schwinger A\ ' = 1/2 [3]. The mass- 
dependent coefficients A 2 and A3 are of higher order; 
the goal of this letter is to provide precise numeri- 
cal values for their 0(a 2 ) and 0(cv 3 ) terms. The rel- 
evance of the results and the improvements with re- 
spect to earlier ones will be discussed separately for 
each lepton. All results were derived using the lat- 
est CODATA [4] recommended mass ratios: m e /m^ = 



4.836 33167(13) x 1CT 3 , m e /m T = 2.875 64(47) x 1(T 4 , 
to m /to g = 206.768 2838(54), m^/m T = 5.945 92(97) x 
1(T 2 , m T /m e = 3477.48(57), to t /to m = 16.8183(27). 
(The value for m T adopted by CODATA in ref. [4] (to t = 
1776.99(29) MeV) is based on the 2002 PDG result [5]. 
This PDG result remains unchanged to date [6, 7].) 



II. ELECTRON 



A. Two-loop contributions 



Seven diagrams contribute to the fourth-order coeffi- 
cient Af \ one to A^ {m e /m^) and one to A^ (m e /m T ). 
As there are no two-loop diagrams contributing to 
a QED £ na t contain both virtual muons and taus, 

A 3 ^ ' (m e /m^,m e /m T )~Q. The mass-independent coeffi- 
cient has been known for almost fifty years [8]: 



A\ 



(4) 



197 3 _ tt 

144 + 12 + 4^ 2 
= -0.328 478 965 579193 78 
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where £(s) is the Riemann zeta function of argument s. 
The coefficient of the two-loop mass-dependent contribu- 
tion to af ED , A 2 4 \l/x), with x — nij/mi, is generated 
by the diagram in fig. 1, where j is the virtual lepton in 
the vacuum polarization subgraph. This coefficient was 




FIG. 1: The QED diagram generating the mass-dependent 
part of af ED in order a 2 . 



first computed in the late 1950s for the muon g— 2 with 



2 



x = m e /m^ -C 1, neglecting terms of 0{x) [9]. The ex- 
act expression for < x < 1 was reported by Elend in 
1966 [10]. However, its numerical evaluation was con- 
sidered tricky because of large cancellations and difficul- 
ties in the estimate of the accuracy of the results, so 
that common practice was to use series expansions in- 
stead [11-13]. Taking advantage of the properties of the 
dilogarithm Li 2 (z) = — f*(dt/t) ln(l — t) [14], the exact 
result was cast in [15] in a very simple and compact ana- 
lytic form, valid, contrary to the one in [10], also for x > 1 
(the case relevant to a® ED and part of a^ ED ): 

Af\l/x) = - j| - ^ + x 2 (4 + 31H*) +|(1 - 5* 2 ) x 



7T 2 fl-X 

— In a; In I - Li 2 (x) + Li 2 (-a;) 

2 \ 1 + x ' 

- — - 2\nx\n(- - x) -Li 2 (x 2 ) 
3 V x 



+ 
(4) 



Forz = l, cq. (4) gives ^ 4) (1) = 119/36-tt 2 /3; of course, 
this contribution is already part of in cq. (3). Nu- 
merical evaluation of eq. (4) with the mass ratios given 
in sec. I yields 



4 4 VeH) 

A^\m e /m T ) 



5.197 386 70(28) x 10~ 7 (5) 
1.837 62 (60) x 10" 9 , (6) 



where the standard errors are only due to the uncertain- 
ties of the mass ratios. The results of eqs. (5) and (6) are 
equal to those obtained with a series expansion in powers 
of y and Iny, with y<^l [4]. 

Adding up eqs. (3), (5) and (6) we get the two-loop 
QED coefficient 



= 4 4) +4 4 Ve/"V)+4Ve/m T ) 

= -0.328 478 444 002 90(60). 



(7) 



The mass-dependent part of Ci^ is small but not negligi- 
ble, giving a relative contribution to the theoretical pre- 
diction of the electron g— 2 of 2.4 parts per billion (ppb). 
This value is much larger than the fabulous 0.7 ppb rel- 
ative uncertainty very recently achieved in the measure- 
ment of a e [2]. The uncertainties in if'^e/m,,) and 
4^ (m e /m T ) are dominated by the latter and were added 
in quadrature. The resulting error 6C { e 4) = 6xl0" 13 leads 
to a totally negligible O(10~ 18 ) uncertainty in the a® ED 
prediction. 



B. Three-loop contributions 



More than one hundred diagrams are involved in the eval- 
uation of the three-loop (sixth-order) QED contribution. 
Their analytic computation required approximately three 
decades, ending in the late 1990s. The coefficient A^ 



arises from 72 diagrams. Its exact expression, mainly 
due to Remiddi and his collaborators, reads [16, 17]: 



A\' = —-n C(3) - — C(5) - -t-^ttt 



+ - 
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810 
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Li 4 (l/2) + — (In 2 2 - 7T 2 ) In 2 2 



1.181241456 587.... 



(8) 



This value is in very good agreement with previous results 
obtained with numerical methods [18]. 

The calculation of the exact expression for the co- 
efficient A^\mi/m,j) for arbitrary values of the mass 
ratio mi/rrij was completed in 1993 by Laporta and 
Remiddi [19, 20] (earlier works include refs. [21]). Let us 
focus on a^ ED (l = e, j = /j,,r). This coefficient can be fur- 
ther split into two parts: the first one, A^ (mi /mj , vac) , 
receives contributions from 36 diagrams containing cither 
muon or r vacuum polarization loops [19], whereas the 

second one, A^\mi/mj, lbl), is due to 12 light-by-light 
scattering diagrams with either muon or r loops [20] . The 
exact expressions for these coefficients are rather compli- 
cated, containing hundreds of polylogarithmic functions 
up to fifth degree (for the light-by-light diagrams) and 
complex arguments (for the vacuum polarization ones). 
Indeed, they were too long to be listed in [19, 20] (but 
were kindly provided by their authors), although series 
expansions were given for the cases of physical relevance. 
The exact expressions for the light-by-light contributions 
also contain a few pentalogarithms in integral form. We 
expressed these integrals in terms of harmonic polyloga- 
rithms (introduced by Remiddi and Vermaseren in [22]), 
thus avoiding their numerical integration. 

The numerical evaluations of the exact expressions for 
4 (mi jrrij , vac) and 4 (mi /rrij , lbl) require some care, 
as the presence of large cancellations makes them prone 
to potentially large roundoff errors. For this reason, nu- 
merical evaluations were carried out with Mathematica 
codes employing exclusively arbitrary-precision numbers, 
keeping track of precision at all points [23]. Harmonic 
polylogarithms were implemented via the Mathematica 
package HPL [24]. Using the recommended mass ratios 
given in sec. I, we obtain the following values: 



j4^(m e /m M ,vac) = 

4 6) (m e /m„,lbl) = 

A^\m e /m T , vac) = 

4 6) (TO e /m r ,lbl) = 



-2.176 84015 (11) x 10~ 5 (9) 
1.439 445 989(77) x 10~ 5 (10) 

-1.167 23(36) x 10" 7 (11) 
5.0905(17) x 10" 8 . (12) 



The sums of eqs. (9)-(10) and eqs. (11)-(12) are 



'(TOe/TO^) 

4 6 '(m e /m T ) 



-7.373 94164(29) x 10~ 6 (13) 
-6.5819(19) x 10" 8 . (14) 
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Equations (9)-(14) provide the first evaluation of the full 
analytic expressions for these coefficients with the CO- 
DATA mass ratios of [4]; they are almost identical to 
the results Af ] (m e /mj = -7.373 941 58(28) x 1(T 6 and 
Af ] {m e /m T ) = -6. 5819(19)xlCT 8 obtained in [4] via the 
approximate series expansions in the mass ratios. The 
small difference between A^\m e /m^) of [4] and eq. (13) 
mainly origins from the 0((m e /m M ) 6 ) term in the series 
expansion of A 2 (me/m^lbl); indeed, due to its small- 
ness, this term was neglected in the expansions [20] used 
in [4]. Expanding the exact Laporta-Remiddi expres- 
sion for the sum of light-by-light and vacuum polarization 
contributions, for r = m(/m 3 <l, we get 



4 6) (r) 

Mr) 
Mr) 

Mr) 

Mr) 



I> 2< 
i=l 



f 2l (r)+0(r w \n 2 r) 

231nr 3((3) 
135 + 2 
4337 In 2 r 



(15) 



He! 
45 



74957 

97200' 



(16) 



22680 476280 
1919tt 2 451205689 



209891 lnr 1811C(3) 



68040 533433600' 
2807 ln 2 r 665641 lnr 
21600 + 2976750 
169 6 77T 2 246800849221 



2304 



3077C(3) 
5760 



(17) 



+ 



907200 480090240000' 



(18) 



55163 ln 2 r 24063509989 lnr 9289((3) 



594000 
340019vr 2 



172889640000 
896194260575549 



23040 



24948000 2396250410400000 



(19) 



The functions M r ) and fi(r) coincide with the expan- 
sions provided in [20], and fe(r) agrees with the com- 
bination of parts from [19] (for the vacuum polariza- 
tion contribution) and [25] (heavy-mass expansions for 
the light-by- light diagrams); M r ) is new - The value of 
{m e /m^) obtained with eq. (15) perfectly agrees with 
that in eq. (13) determined with the exact formulae. In- 
deed, their difference is of O(10~ 23 ), to be compared with 

the O(10~ 13 ) error 8A^(m e /mf^) due to the present un- 
certainty of the ratio me/m^. Therefore, it will be possi- 
ble to compute A^ ^(m e /m M ) with the simple expansion 
in eq. (15) - thus avoiding the complexities of the exact 
expressions - even if the precision of the ratio rrie/m^ 
will improve in the future by orders of magnitude. 

The contribution of the three-loop diagrams with both 
\x and r loop insertions in the photon propagator can be 
calculated numerically from the integral expressions of 
ref. [11]. We get 



4 6) (TO e /m M ,m e /TO T ) = 1.909 45 (62) x 10 



-13 



(20) 



a totally negligible O(10 21 ) contribution to a|? ED . 
Adding up eqs. (8), (13), (14) and (20) we obtain the 



three-loop QED coefficient 

Cf ] = 1.181234016 827(19). 



(21) 



The relative contribution to a® ED of the mass-dependent 
part of this coefficient is ~ 0.1 ppb. This is smaller than 
the present ~ 0.7 ppb experimental uncertainty [2]. The 
error 1.9 xlO -11 in eq. (21) leads to a totally negligible 
O(10" 19 ) uncertainty in a« ED . 



C. Determination of a. from the electron g — 2 

Very recently, a new measurement of the electron anoma- 
lous magnetic moment by Gabrielse and his collaborators 
achieved the fabulous relative uncertainty of 0.66 ppb [2], 



a EXP = 115 965 218 0.85 (76) x 10 



-12 



(22) 



This uncertainty is nearly six times smaller than that 
of the last measurement of a e reported back in 1987, 
a EXP = 1159652188.3(4.2)xl0- 12 [4, 26]. These two mea- 
surements differ by 1.7 standard deviations. 

The fine-structure constant a can be determined 
equating the theoretical SM prediction of the electron 
q— 2 with its measured value 



(23) 



The SM prediction contains the QED contribution 
a^ ED (a) = J2i=i C'e 2l \a/TTY (higher-order coefficients 
are assumed to be negligible), plus small weak and 
hadronic loop effects: a b e M (a) = a« ED (a)+a EW +a« AD (the 
dependence on a of any contribution other than a® ED is 
negligible). The electroweak contribution is [4]: 



a EW = 0.0297(5) x 10 



-12. 



(24) 



this precise value includes the two-loop contributions cal- 
culated in ref. [27]. The hadronic term is [4, 28]: 



< AD = 1.671 (19) x 10 



-12 



(25) 



The latest value for the four-loop QED coefficient is 
C*i 8) = -1.7283(35) [29] . Following the argument of [4] , we 
adopt the educated guess Ci 10 ^= 0.0(3.8) for the five-loop 
coefficient. The errors 8C ( e S) = 0.0035 and SC ( e W) = 3.8 
lead to an uncertainty of 0.1 x 10~ 12 and 0.3 x 10~ 12 in 
aJ? ED , respectively. Solving eq. (23) with the new mea- 
sured value of eq. (22), we obtain 



= 137.035 999 709(12) (30) (2) (90) 
= 137.035 999 709 (96) [0.70 ppb]. 



(26) 



The first and second errors are due to the uncertainties of 
the four- and five- loop QED coefficient 5C^ and <5Ce 10 \ 
respectively; the third one is caused by the tiny i5ag AD , 
and the last one (90 x 10~ 9 ) is from the experimental 
<5a EXP in eq. (22). The uncertainty of the electroweak and 
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two/three-loop QED contributions are totally negligible 
at present. The determination in eq. (26) is in perfect 
agreement with the new result of ref. [30] , 



a' 1 = 137.035 999 710(96) 



(27) 



(also based on the new measurement of ref. [2]), whose 
great precision represents the first significant improve- 
ment of this fundamental constant in a decade. The to- 
tally negligible difference between eqs. (26) and (27) is 
due to the rounded value a EW = 0.030(1) x 10~ 12 [4] 
employed by the authors of ref. [30] instead of eq. (24). 

At present, the best determinations of a independent 
of the electron g— 2 are 

a _1 (Rb) = 137.035 998 78 (91) [6.7 ppb], (28) 
aT 1 (Cs) = 137.036 000 0(11) [8.0 ppb]; (29) 

they are less precise by roughly a factor of ten. The value 
a _1 (Rb) was deduced from the measurement of the ra- 
tio h/Mnb based on Bloch oscillations of Rb atoms in 
an optical lattice (h is the Planck constant and M^b is 
the mass of the Rb atom) [31], while a _1 (Cs) was deter- 
mined from the measurement of the ratio h/Mc s (Mq s 
is the mass of the Cs atom) via cesium recoil measure- 
ment techniques [32, 33]. These two determinations of 
a also rely on the precisely known Rydberg constant 
and relative atomic masses of the electron, Rb and Cs 
atoms [4, 34]. The values of a in eqs. (28) and (29) are 
in good agreement with the result of eq. (26), differing 
from the latter by —1.0 and +0.3 standard deviations, 
respectively. This comparison provides a beautiful test 
of the validity of QED. It also probes for possible electron 
substructure [30]. 



III. TAU 

The two-loop mass-dependent QED contributions to the 
anomalous magnetic moment of the r, obtained by direct 
evaluation of the exact formula in eq. (4), are 



A^\m T /m e ) 



2.024 284 (55), 
0.361652 (38). 



(30) 
(31) 



These two values are very similar to those computed via a 
dispersive integral in ref. [35] (which, however, contain no 
estimates of the uncertainties). Equations (30) and (31) 
are also in agreement (but more accurate) with those of 
ref. [12]. Adding up eqs. (3), (30) and (31) we get 



Cj 4) = 2.057457(93) 



(32) 



(note that the uncertainties in m T /m e and m T /m^ are 
correlated). The resulting error 9.3 x 10~ 5 leads to a 
5x 10~ 10 uncertainty in a^ ED . 

We computed the three-loop mass-dependent contribu- 
tions by direct numerical evaluation of the exact analytic 



expressions (see sec. II B). The results are: 

4 6) (m T /m e ,vac) = 7.256 99(41) (33) 

4 6) (m r /m e ,lbl) = 39.1351 (11) (34) 

4 6) (m T /m^,vac) = -0.023 554(51) (35) 

4 6) (m T /m M ,lbl) = 7.033 76 (71). (36) 

Employing the approximate series expansions (sec 
sec. II B) we obtain almost identical values: 7.25699(41), 
39.1351(11), -0.023564(51), 7.03375(71). The estimates 
of ref. [35] were: 10.0002, 39.5217, 2.9340, and 4.4412 (no 
error estimates were provided), respectively; they are at 
variance with our results, eqs. (33)-(36), derived from 
the exact analytic expressions. The estimates of ref. [36] 
compare slightly better: 7.2670, 39.6, -0.1222, 4.47 (no 

errors provided). In the specific case 
it's easy to check that the values of refs. [35, 36] differ 
from eq. (36) because their derivations did not include 
terms of 0(m^/m T ), which turn out to be unexpectedly 
large. The sums of eqs. (33)-(34) and (35)-(36) are 

4 6) (m r /m e ) = 46.3921 (15), (37) 
4 6) KK) = 7.01021 (76). (38) 

The contribution of the three-loop diagrams with both 
electron- and muon-loop insertions in the photon prop- 
agator can be calculated numerically from the integral 
expressions of [11]. We get 

Af l {m T /m e ,m T /m lJ ) = 3.34797(41). (39) 

This value disagrees with the results of refs. [36] (1.679) 
and [35] (2.75316). Combining the three-loop results of 
eqs. (8), (37), (38) and (39) we find the sixth-order QED 
coefficient 



Cf> = 57.9315(27). 



(40) 



The error 2.7 x 10~ 3 induces a 3 x 10~ n uncertainty in 
fl QED or( j er Q f ma g n itude of the three-loop contribu- 
tion to a^ ED , dominated by the mass-dependent terms, is 
comparable to that of electroweak and hadronic effects. 

Adding up all the above contributions and using the 
new value a" 1 = 137.035 999 710(96) [30] (or the value 
derived in eq. (26), o^ 1 = 137.035 999 709(96) - the 
difference is totally negligible) we obtain the total QED 
contribution to the g— 2 of the r: 



117324 (2) x 10" 



(41) 



The error (5a^ ED is the uncertainty SC^P (a/7r) 4 <~ 
7T 2 ln 2 (m T /m e )(a/7r) 4 ~ 2 x 10~ 8 which we assigned to 
a QED f or unca i cu i a ted four- loop contributions. As we 
mentioned earlier, the errors due to the uncertainties of 
the 0(a 2 ) and 0(a 3 ) terms are negligible. The error in- 
duced by the uncertainty of a is only 8xl0~ 13 (and thus 
totally negligible). 

The g— 2 of the r is a very interesting observable, even 
if the short lifetime of this lepton makes its measurement 
very difficult at present. The possibility to improve the 
recent experimental bounds [7] is certainly not excluded. 
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IV. MUON 

This final section reports the results relevant to a^ ED (see 
[37] for recent reviews of the entire SM prediction). Some 
of them were already presented in [15]. The two-loop 
contributions are 

A^im^/me) = 1.094 258 3111 (84), (42) 
A^\m^/m T ) = 0.000 078 064(25). (43) 

The sum of eqs. (3), (42) and (43) provides the coefficient 
Cf =0.765 857 410 (27). The value SCjp = 2.7xl0~ 8 was 
obtained adding in quadrature the errors in eqs. (42) and 
(43). It produces a tiny 1.4 xlO -13 uncertainty in a^ ED . 



The three-loop contributions 


are 




^2 6) ("V/ TO e: vaC ) = 


1.920 455130 (33) 


(44) 


Af\mJm e M) = 


20.947 924 89(16) 


(45) 


A^ (m^/m r , vac) = 


-0.001782 33 (48) 


(46) 


Af\mJm T M) = 


0.002142 83(69) 


(47) 


Af\m^/m e ) = 


22.868 380 02 (20) 


(48) 


Af\m^/m T ) = 


0.000 360 51 (21). 


(49) 



The analytic calculation of the three-loop diagrams with 
both electron and r loop insertions in the photon prop- 
agator became available in 1999 [13] and was confirmed 
more recently [38] . This analytic result yields the numer- 
ical value [15] 



Af'im^/m^m^/mr) = 0.000 52766(17), (50) 

providing a small 0.7 x 10~ n contribution to a^ ED . The 
error 1.7 x 10~ 7 is caused by the uncertainty of the ratio 
"V/ m r- Combining the three- loop results of eqs. (8), 
(48), (49) and (50) we get the three-loop coefficient 
C^ 6) = 24.050 509 64(43). The error 4.3 x 10~ 7 induces 
a negligible O(10~ 14 ) uncertainty in a^ ED . 

Adding the four-loop and leading five-loop contri- 

butions computed by Kinoshita and Nio, Cy, = 

130.9916 (80) [29, 39] and C^ 10) = 663 (20) [40] (estimates 
obtained with the renormalization-group method agree 
with this five- loop result [41]), and using the new value 
a- 1 = 137.035 999 710(96) [30] (or the value derived in 
eq. (26), a- 1 = 137.035 999 709(96) - the difference is 
negligible) we get the new total QED contribution to the 
muon g— 2, 

a« ED = 116 584 718.09 (14) (08) x HT 11 . (51) 



The first error is determined by the uncertainties of 
the QED coefficients (dominated by the five-loop one, 

— 20), while the second is caused by the tiny 
uncertainty 8a. Equation (51) is in good agreement 
with the recent value a« ED = 116 584 717.62 (14) (85) x 
10~ n [40], and the uncertainty due to 5a is strongly re- 
duced. 

Note added after publication 

The value of the mass-independent eighth-order QED 
coefficient A\ has been recently revised by Kinoshita 
and collaborators [42], inducing the following updates. 

The revision of A± (and, consequently, of the 
total four-loop QED coefficient — A^), from 

-1.7283(35) [29] to -1.9144(35) [42], induces a shift in 
the value of a from eq. (26) to 

a' 1 = 137.035 999 068(12) (30) (2) (90) 

= 137.035 999 068(96) [0.70 ppb]. (52) 

This new value is still in good agreement with those in 
eqs. (28) and (29), which are less precise by roughly a 
factor of ten and differ from the value in eq. (52) by —0.3 
and +0.8 standard deviations, respectively. 

The total QED contribution to the muon g— 2 shifts 
from the value in eq. (51) to 

a« ED = 116 584 718.10 (14) (08) x HT 11 . (53) 

This tiny variation is only of 0((a/7r) 5 ). Indeed, the 
0((a/7r) 4 ) shift in a« ED due to the revision of Af ] (and, 
consequently, of the total four-loop QED coefficient, now 
standing at C/? = 130.8055 (80)) is compensated by the 
0{{a/-K) 4 ) change in the value of a determined from the 
electron g— 2. 

The shift in the value of a from eq. (26) to eq. (52) 
induces no appreciable variation in the total QED contri- 
bution to the g— 2 of the r lepton. 
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